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Magneto-optical spectroscopy of (Ga,Mn)N epilayers 
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We report on the magneto-optical spectroscopy and cathodoluminescence of a set of wurtzite 
(Ga,Mn)N epilayers with a low Mn content, grown by molecular beam epitaxy. The sharpness 
of the absorption lines associated to the Mn 3+ internal transitions allows a precise study of its 
Zeeman effect in both Faraday and Voigt configurations. We obtain a good agreement if we assume 
a dynamical Jahn- Teller effect in the 3d 4 configuration of Mn, and we determine the parameters 
of the effective Hamiltonians describing the 5 T2 and E levels, and those of the spin Hamiltonian 
in the ground spin multiplet, from which the magnetization of the isolated ion can be calculated. 
On layers grown on transparent substrates, transmission close to the band gap, and the associated 
magnetic circular dichroism, reveal the presence of the giant Zeeman effect resulting from exchange 
interactions between the Mn 3+ ions and the carriers. The spin-hole interaction is found to be 
ferromagnetic. 

PACS numbers: 75.50.Pp, 75.30.Hx, 78.20.Ls 
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I. INTRODUCTION 

Extrapolating the Zener model of carrier-induced fer- 
romagnetism to wide bandgap diluted magnetic semicon- 
ductors (DMS) predicts that high critical temperatures 
should be achieved, provided several demanding assump- 
tions are satisfied. In particular, p-type (Ga,Mn)N would 
be ferromagnetic with a critical temperature exceeding 
room temperature. This implies however the incorpora- 
tion of 5% of Mn into GaN, substituting Ga in the form 
of Mn ions, and a strong p-type doping. 1 

Early experimental works following this prediction did 
not give any clearcut conclusion about the magnetic 
properties of (Ga,Mn)N: ferromagnetic properties at high 
temperature^ as well as paramagnetic properties down 
to low temperature^^ have been reported. However, 
there is now a general agreement that (Ga,Mn)N does 
not offer the configuration which exists in (Ga,Mn)As, 
with the transfer of an electron from the valence band to 
the d- levels, resulting in both the d 5 configuration (with 
a spin 5/2) which is observed and the acceptor character 
of the Mn impurity. 

Ab initio calculations of the band structure 6 ' 7 ' 8 ' 9 ' 10,11 
do not exclude the possibility of a ferromagnetic interac- 
tion, but they all conclude that the d bands of Mn are 
located rather high in the bandgap of GaN, although its 
exact position cannot be given. 

From an experimental point of view, diverging conclu- 
sions have been drawn about the valence state of Mn 
in (Ga,Mn)N. In bulk (Ga,Mn)N with a very low Mn 
content, the d 5 state was detected by electron paramag- 
netic resonance^ in similar samples, the d 4 valence state 
was deduced from magneto-optical measurements upon 



co-doping with Mg£ The Mn content in these samples 
was chosen to be very low, a few 10 18 cm~ 3 at most, 
so that these two results can be understood from the 
strong n-type character of the first samples and the usual 
acceptor character of Mg in GaN. The same d 5 state 
was observed by x-ray absorption spectroscopy at the 
L-edge of Mn near the surface of layers grown by plasma- 
assisted molecular beam epitaxy on n-type, Sn-doped 
GaN templates, 12 while we have observed the d A state 
by x-ray absorption spectroscopy at the K-edge of Mn in 
layers grown also by plasma-assisted molecular beam epi- 
taxy but on undoped GaN templatespii^i 3 - Various types 
of valence states, from d 5 to d 3 , were invoked to explain 
optical spectra of epilayers with a larger Mn contentii 4 ^ 
These results were mainly obtained on (Ga,Mn)N with 
the wurtzite structure : the case of cubic (Ga,Mn)N lay- 
ers appears to be peculiar, since p-type doping due to the 
incorporation of Manganese was reportedfi 6 ^ 7 -. 

Optical absorption spectr gi^l&lS^ show (i) a struc- 
tured absorption band in the near infra-red, with a sharp 
zero-phonon line, which was attributed to the d—d transi- 
tions of the Mn 3+ ion (the d A configuration), (ii) a broad 
absorption band attributed to the transitions between 
the Mn levels and the bands of the semiconductor, (in) 
absorption at the bandgap of the semiconductor. 

The present work is devoted to a magneto- 
spectroscopic study of (Ga,Mn)N samples with a low Mn 
content, so that quite sharp features are observed, and to 
the evolution of the magnetic and magneto-spectroscopic 
properties upon increasing the Mn content. Main re- 
sults are a determination of the parameters which enter 
the effective Hamiltonian of the ground level of Mn in 
(Ga,Mn)N, which results from a dynamical Jahn- Teller 



effect in the d 4 configuration of Mn, and the observation 
of spin-carrier coupling through the giant Zeeman effect 
at the bandgap of (Ga,Mn)N. 



II. SAMPLES AND EXPERIMENTAL SET-UP 

Epilayers of wurzite (Ga,Mn)N were grown by 
nitrogen-plasma-assisted molecular beam epitaxy 
(MBE), on a GaN template, previously grown by metal- 
organic chemical vapor deposition (MOCVD) on the c 
surface of a sapphire substrate. A few additional samples 
were grown on A1N templates, so that the substrate is 
transparent in the vicinity of the bandgap of (Ga,Mn)N. 
Unfortunately, the growth had to be initiated with a 
GaN layer, which was kept as thin as possible. The 
growth temperature was kept at 720° C, which is our 
usual growth temperature for GaN. A special attention 
was paid to the different growth regimes related to the 
Ga/N flux ratio and to the manganese fluxi^i Samples 
were characterized in situ using reflection high-energy 
electron diffraction (RHEED). Post-growth characteriza- 
tion comprised ex situ secondary ion mass spectroscopy 
(SIMS) and x-ray diffraction (XRD). 

Optical absorption in the near infrared was measured 
using a Fourier Transform Infrared spectrometer (FTIR) 
with a spectral resolution of 0.5 meV, using a tungsten 
lamp as the light source and a silicon photodetector. 
Magneto-optical spectroscopy in the near infrared and in 
the vicinity of the band gap energy was measured using 
a CCD camera attached to a grating spectrometer. The 
magnetic circular dichroism was then calculated from the 
transmission spectra of right- and left-circular polarized 
light. In some cases, it was also measured directly using 
a 50 kHz photoelastic modulator and a photomultiplier 
tube. 

Cathodoluminescence was observed in a scanning elec- 
tron microscope, with the sample placed on a cold finger 
at 5 K, using an accelerating voltage of 4 kV and a typical 
current of the order of 30 nA. The luminescence intensity 
was recorded on a GaAs photomultiplier tube attached 
to a 460 mm spectrometer for imaging, and a charge cou- 
pled device (CCD) camera for spectroscopy. Attempts to 
measure the photoluminescence excited by a He-Cd laser 
(20 mW on a 50 /^m spot size) were unsuccessful. 

All the transmission spectra exhibited interference 
fringes due to internal reflexions between the surface and 
the sapphire/GaN interface. These interferences were re- 
moved using a lineshape fit of the transmittance which 
includes a sin(2ne/A) function, where e is the thickness of 
the GaN + (Ga,Mn)N layer, A the wavelength, and n the 
optical index of GaN at room temperature^ This pro- 
cess is quite efficient in the case of sharp lines, such as the 
near infrared zero-phonon line, for which the fit can be 
optimized over a narrow range across the line. The pro- 
cess is less accurate for broad absorption bands (broader 
than 25 meV typically), particularly in the vicinity of 
the bandgap, since we should use the optical index of 



(Ga,Mn)N at low temperature. 

Magnetization measurements were performed using a 
superconducting quantum interference device (SQUID) 
magnetometer at magnetic fields up to 5 T applied in 
the plane of the sample (_L c axis) or perpendicular to 
the plane (|| c axis), at temperatures down to 2 K. The 
large diamagnetic contribution from the substrate was 
evaluated on a piece of the same substrate. 



III. MAGNETO-OPTICAL SPECTROSCOPY OF 
d-d TRANSITIONS 

In this section, we present the spectroscopic study 
in the near infrared, including magneto-transmission 
and cathodoluminescence. Previous transmission spec- 
tra have been attributed^ to the 5 T2 — ► 5 E transition of 
the d A configuration of Mn, involving a large cubic crys- 
tal field splitting, about 1.4 eV. The better resolution 
achieved on the present samples with a low Mn content, 
and the comparison between transmission and cathodo- 
luminescence spectra, lead us to reconsider the effective 
Hamiltonian which takes into account the Jahn- Teller ef- 
fect and the trigonal crystal field. 



A. Experimental results 

Fig. n displays a transmission spectrum and a 
cathodoluminescence spectrum of a (Ga,Mn)N layer with 
0.06% Mn. A series of cathodoluminescence spectra were 
measured locally from different points of the sample. 
The one displayed in Fig. ^ is typical of what was ob- 
tained from most points ; occasionally, broadened spec- 
tra were observed. The origin of the horizontal scale in 
Fig. His the common position of the zero-phonon line at 
1413 meV, and the scale is inverted for the two spectra. 

Most features are present on both spectra : they are 
clearly phonon replica. They may involve local modes 
of the Mn impurity and its cluster of nearest neigh- 
bors, which are expected around 20 meV and 70 meV. 18 
These frequencies are close to high symmetry phonons 
of GaNj2ii24i2ii as expected since the difference of mass 
between Mn and Ga is small. The presence of resonances 
with the phonons of GaN makes the spectrum quite com- 
plex, but the coupling to the mode at 20 meV dominates, 
so that a simple model of coupling to a single local mode 
can be a reasonable approximation. 

The coupling to this principal mode is rather weak: 
the intensity of the zero-phonon line and that of the 
phonon replica are of the same order of magnitude. In 
the case of a single local mode, the intensity of the n th 
replica to the zero-phonon line is expected to be S n /n\, 
where S is the so-called Huang-Rhys factor: from the first 
replica, we obtain S =0.7 from the transmission spec- 
trum, S =1.9 from the cathodoluminescence spectrum. 
This is in good agreement with previous determinations 
(S =0.6 in Ref.Gj, S =1.1 in Ref-EK). The main contri- 
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FIG. 1: Absorption and cathodoluminescence spectra of a 
(Ga,Mn)N layer with 0.06% Mn. The horizontal scale is the 
distance from the zero phonon line at 1413 meV. The posi- 
tion of high-symmetry phonons of GaN— a— <— is indicated. 
Local modes around the Mn impurity are expected at 20 and 
70 meV. Diamond and star indicate additional zero-phonon 
lines. 



bution to the line at AE—40 meV arises from the second 
phonon line due to the same mode. 

This coupling is much weaker than in the case of Cr 
in II- VI semiconductors r& which has the same d 4 con- 
figuration. Then the zero-phonon line is weak (3.7 x 
10 -27 m 2 eV/Cr in ZnSe^ to be compared to 1.1 x 
1CT 23 m 2 eV/Mn in GaN). The intensity is transferred 
to a broad vibronic band, three orders of magnitude more 
intense than the zero-phonon line (so that the total in- 
tegrated intensity is about the same in both cases). The 
Huang-Rhys factor, as deduced from the position of the 
maximum of the vibronic band, at SHlj from the zero- 
phonon line, is as high as 13 for Cr 2+ in ZnTe and 18 in 
CdS. 

Two lines in Fig. ^ cannot be ascribed to phonon 
replica. The line marked with a diamond in Fig. ^ is 
present only in the cathodoluminescence spectrum. It 
must correspond to a transition to an excited state of 
the ground 5 T2 multiplet. Another line, marked with 
a star in Fig. ^ is observed in transmission but not in 
cathodoluminescence. Its intensity is much lower than 
that of the zero-phonon line at 1413 meV, and we ten- 
tatively ascribe it to a spin-forbidden transition to one 
of the 3 Ti or 1 Ai levels, which are indeed quite close for 
such a value of the cubic crystal fieldiS 7 . 

Fig. [3 shows a series of transmission spectra in the 
range of the zero-phonon line, for another sample with a 
low Mn content (0.03%), for different values of field and 
temperature, in both the Faraday and Voigt configura- 
tions. The general features reported by Wolos et ali^ are 
observed, with however a higher resolution. Spectra up 
to 22 T have been also recorded from the same sample, 
with a slightly lower resolution (not shown). 

Two points should be stressed. First, the zero-field 
spectrum contains essentially one line : any structure, 



if present, should be comprised within the observed 
linewidth of 1.1 meV. Second, the "shift" at 7 T reported 
in Ref. 14 is actually an intensity transfer between two 
lines labeled LI and L2. From a plot of the positions and 
intensities deduced from a Gaussian fit (Fig. |3J) , a third 
peak L3 appears at intermediate magnetic field, which 
exhibits a small but systematic shift with respect to L2. 

The appearance of L3 at 3 T, on one hand, and the 
intensity transfer to L2 observed between 10 and 13 T, 
on the other hand, indicate two values of the applied field 
where a crossing takes place in the manifold of ground 
states. The relative position of the ground states which 
serve as initial states for the absorption lines can be de- 
duced from spectra recorded as a function of the temper- 
ature for selected values of the magnetic field (Fig. 0}. 
For instance, it is clear from the spectra that two lines, 
at 1411.8 meV and at 1412.8 meV, co-exist even at low 
temperature at 6 T. At lower fields only one line is seen 
at 1.7 K and a second line rises on its low energy side at 
higher temperature ; the opposite is observed at higher 
fields. A plot of the intensities (not shown) is well re- 
produced by an activation energy with Lande factors of 
the order of 2. As the positions of LI, L2 and L3 remain 
almost unchanged (at most 0.2 meV over a 10 T scan), 
almost identical Zeeman shifts occur in both the ground 
and excited multiplets - which may be accidental, but 
rather suggests that the spin-orbit coupling is weak. 

This is supported by the fact that we measure a weak 
(although non zero) magnetic circular dichroism (MCD), 
[I(a + )-I(a-)]/[I(a + )+I{a-)]. Fig. \Sjp displays spectra 
recorded for the two circular polarizations, er + and <r~, 
at 11 T and 1.7 K, together with the resulting MCD 
(Fig. |SJi) . Measuring such a small dichroism on a thin 
layer is difficult. MCD spectra with a much higher signal- 
to-noise ratio, obtained on bulk material, are given in 
Ref. UjJ. The same lineshape is observed at high field, 
but the MCD persists at lower field in Ref. [lj while we 
observe a faster decrease and even a trend to a change of 
sign (Fig. Ob) . 

We shall show now, by constructing an effective Hamil- 
tonian, that these spectra are well explained by a d 4 con- 
figuration, in a crystal field of trigonal symmetry, with a 
dynamical Jahn- Teller effect which reduces the spin-orbit 
coupling. 



B. Effective Hamiltonian 

The Hamiltonian of a d 4 ion contains the crystal field 
terms (cubic and trigonal components), the spin-orbit 
coupling, the Zeeman effect, the coupling to local strain, 
and the Jahn- Teller effect. 

The cubic component of the crystal field splits the 5 D 
multiplet into a ground 5 T2 orbital triplet and a 5 E or- 
bital doublet (noted lODq in table I). A basis of orbital 
states of the triplet adapted to the cubic symmetry is 
given by the three states |£ >, \r) >, |? > which trans- 
form like yz, zx, xy where x, y, z label the cubic axes; 
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FIG. 2: (color online) Transmission spectra of a (Ga,Mn)N layer with 0.03% Mn, for different fields applied in the Faraday 
configuration at 1.7 K (a) and 10 K (b), and in the Voigt configuration at 4.2 K (c). Black lines show the experimental spectra 
and red lines the calculated ones. 



for the doublet, a basis is given by \8 >, \e >. 

In the case of a wurtzite crystal, there is also a trigonal 
crystal field. It has off-diagonal matrix elements in the 
cubic basis of the orbital triplet and no matrix elements 
within the doublet. In the case of epitaxial layers grown 
along the c axis, it may contain a contribution from the 
biaxial strain due to the partially compensated lattice 
mismatch. 

For the free Mn 3+ ion, the spin-orbit coupling within 
the 5 D multiplet is written AL.S with L = 2, S = 2 and 
A w 11 meV 3ln a crystal, this value can be reduced 
slightly by hybridization, and significantly by the Jahn- 
Teller effect (see below). 

The coupling to local strain (" strain fluctuations" ) in- 
volves trigonal and tetragonal components. The coupling 
to trigonal strain, i.e., strain components of T2 symme- 
try, has off-diagonal matrix elements in the cubic basis 
of the orbital triplet and no matrix elements within the 
doublet. The coupling to tetragonal strain, i.e., strain 
components of E symmetry, has matrix elements within 
the triplet (where they are diagonal) and the doublet. 

Finally, there is a Jahn- Teller coupling of the orbital 
triplet to E and T2 modes, and of the orbital doublet to 
the E modes. 

The most simple case of Jahn- Teller coupling is that 
of an orbital triplet coupled to a single E mode. A very 



intuitive treatment, based on first and second order per- 
turbations, was given by Ham 29 . Basically, the orbital 
triplet is replaced by a vibronic triplet : each orbital 
state is associated to a potential well with its minimum 
corresponding to a distortion of the environment along a 
cubic direction. The vibronic states are the product of 
the orbital state by the eigenstates of the displaced har- 
monic oscillator. The Hamiltonian acting within the or- 
bital triplet, as defined above, is replaced by an effective 
Hamiltonian acting within the vibronic triplet: the main 
result is that the first-order contribution of off-diagonal 
operators is reduced by the overlap of eigenstates of the 
different harmonic oscillators. First-order contributions 
of diagonal operators are not reduced. Second-order con- 
tributions have to be considered : they are not reduced 
as the first-order ones so that their role can be enhanced. 

In the case of a crystal of cubic symmetry, a strong 
Jahn- Teller effect results in a strong reduction of all con- 
tributions but that of tetragonal strain fluctuations. One 
can use the approximation of the " static Jahn- Teller cou- 
pling" : the local tetragonal strain slightly lowers the en- 
ergy of a well corresponding to a distortion along one 
of the three cubic directions, and for this ion all cal- 
culations are performed with this particular distortion. 
This approximation was used in the first study of Cr 2+ in 
II- VI semiconductors, 26 with the additional assumption 
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FIG. 3: (color online) Position of the absorption lines ob- 
served in the (Ga,Mn)N layer with 0.03% Mn, as a function 
of the magnetic field applied in the Faraday and the Voigt 
configurations. The size of the symbols is proportional to 
the intensity. The sample is the same as in Fig. [5] but the 
experiments were realized in a different set-up. 
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FIG. 5: (color online) (a) magnetic circular dichroism for dif- 
ferent values of the applied field (field step lT)(b) Transmis- 
sion of right- and left-polarized light (green and blue lines 
respectively) and (c) magnetic circular dichroism (black line) 
of the (Ga,Mn)N layer with 0.03% Mn, at 1.7 K at 11 T in the 
Faraday configuration. Red line shows the calculated MCD 
spectrum. 
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FIG. 4: (color online) Transmission spectra of the (Ga,Mn)N 
layer with 0.03% Mn, for different fields applied in the Fara- 
day configuration, and different temperatures, as indicated. 
Black lines show the experimental spectra and red lines the 
calculated ones. 



that there is no Jahn- Teller coupling in the 5 E excited 
state. The same assumption was used later— in the case 
of Cr 2+ in CdS, which has the wurtzite structure. More 
recently^ a reasonable agreement with the available data 
was obtained in the case of Mn 3+ in GaN. 

In the present study, two points however cannot be 
explained by this model: (?) the presence of a second 
zero-phonon line in the cathodolumincsccncc spectrum, 
and (ii) the observation of a main zero-phonon line at 
zero field, both in cathodoluminescence and in transmis- 
sion, with a linewidth of only 1 meV. Spectra calculated 
within the static limit of the Jahn- Teller effect (see Fig. 8 



of Ref. ll4Jl show a doublet which could fit the available 
spectra with a broad linewidth, but is ruled out by the 
spectra of the present study (Fig. [I}. We could not obtain 
a better fit with the same model and different parame- 
ters. This leads us to re-examine the treatment of the 
Jahn- Teller effect in the case of (Ga,Mn)N. 

We may note that an improvement of the available 
spectra for Cr 2+ in ZnS and ZnSe recently lead Bevilac- 
qua et o/— to propose a rather weak Jahn- Teller cou- 
pling to the T 2 modes (and not the E mode) for the or- 
bital triplet, and a stronger coupling to the E mode for 
the orbital doublet (instead of no coupling) . Both cases 
are more difficult to handle, and necessitate a more com- 
plete treatment of the vibronic states 31 which is beyond 
the scope of the present study. 

Actually, we obtained a reasonable, although not per- 
fect, agreement in the frame of a dynamic Jahn- Teller 
effect for both the triplet and the doublet, and not a 
static one for Cr in II-VI's. This is not unexpected for 
two reasons: (i) the smaller value of the Huang- Rhys fac- 
tor suggests that the Ham reduction factor— should be 
smaller (this is but a hint, since the Huang- Rhys factor 
measured on spectra involves a vibrational overlap be- 
tween a distorted harmonic oscillator of the triplet and 
a distorted harmonic oscillator of the doublet, while the 
reduction factors involve an overlap between two differ- 
ent harmonic oscillators of two states of the same mul- 
tiplct); (ii) the trigonal crystal field due to the wurtzite 
structure tends to mix the three cubic states and, even 
if it is reduced by the Jahn- Teller effect, it may remain 
larger than the tetragonal component of strain fluctua- 
tions (which tends to stabilize one cubic distortion). 

The effective hamiltonian which must be diagonalized 
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proportional to the matrix element of the appropriate 
projection of L. Spectra were calculated by adding all 
contributions: 



Magnetic field [T] 

FIG. 6: Energy levels calculated taking into account the cubic 
and trigonal crystal field, the spin-orbit coupling, the Zeeman 
effect in both Faraday and Voigt configuration and the Jahn- 
Teller effect. The arrows indicate the main components of the 
optical transitions. The red dotted lines correspond to a fit of 
the energy levels of the ground state using the effective spin 
Hamiltonian (4). The corresponding fit parameters are given 
in table II 



within the 5 T? vibronic multiplet thus reads: 



H 



KV trl [ZL 2 c - L(L + 1)] - kAL.S + /j, b B.(-kL + 2S) 
+p(L x S x L y Sy +perm.) + d{L 2 x Sl. + perm.) 



(1) 



where Vtri measures the trigonal component of the crys- 
tal field, L c is the projection of the pseudo-kinetic mo- 
ment operator L — 1 on the c axis, k is the Ham reduc- 
tion factor which applies on all first-order contributions 
which are off-diagonal in the cubic basis, and the two 
last terms arise from second-order contributions of spin- 
orbit coupling. We assume that strain fluctuations have 
a negligible effect. 

The case of the Jahn- Teller effect in the doublet is 
more complicated. The full treatment for Cr 2+ in ZnS 
and ZnSe^L concludes that there is a significant Jahn- 
Teller effect, and that the excited levels which contribute 
to the absorption lines are comprised within a narrow 
energy range (less than 0.3 meV). We keep a simple de- 
scription and assume a vibronic doublet, with an effec- 
tive Hamiltonian parameterized by d! which contains a 
second-order contribution of spin-orbit coupling (arising 
from the AtbL-S spin-orbit coupling between the ground 
triplet and the orbital doublet), and the Zeeman effect. 
Using the operators XJg and U e defined by H&mM within 
the orbital doublet, we have 

H = d'[U e (2S 2 z - Si - S 2 ) + U E .V3(S 2 X - S^)} 

+^ B B.2S (2) 

The intensity of the 5 T2 — ► E transitions is propor- 
tional to the square of the matrix element of the elec- 
tric dipole, which in a tetrahedral configuration is itself 



/+ = 



E^/o eM-(E-E if ) 2 /2a 2 ) |(z|L*|/)| 2 n(gQ 

(3) 
For absorption, \i > spans all calculated states of the 
5 T2 multiplet, taking into account the thermal occupancy 
n(Ei)=KT Ei l kBT , and \f > all calculated states of the 
5 E multiplet. The full width at half maximum of the 
Gaussian lineshape (equal to 0.8 meV) and its intensity 
were adjusted at 1.7 K and zero-field, and kept constant 
for all temperature and field values. 

Absorption and emission spectra were calculated using 
the parameters listed in table |U Fig. shows the energy 
levels. The position of the second zero phonon line ob- 
served in cathodoluminescence 6 meV above the main 
line corresponds to the splitting of the orbital triplet by 
the trigonal crystal field. Calculated transmission spec- 
tra are shown by red lines in Fig. [3 0] and [5] All lines 
contain several components: the main contributions are 
shown by arrows in Fig. H3 The observed "jumps" are 
due to the crossing between the spin states of the ground 
orbital singlet. The ^-factors in the ground and excited 
multiplets are close to 2 so that the Zeeman shifts are 
parallel. The position in field of the crossings is directly 
related to the zero-field splitting, due to the combined 
effect of the trigonal crystal field and the spin orbit cou- 
pling. From Fig. [21 and 01 the model of Mn in the d 4 
configuration well describes the field and temperature de- 
pendence of the absorption spectra. The MCD at high 
magnetic field is also well described (Fig.EJ. 

The calculated spectra deviate from our experimental 
findings on two points. First, the intensity of the second 
zero-phonon line in emission (labeled by a diamond on 
Fig.EJ is about twice larger that experimentally observed 
in cathodoluminescence. Second, the calculated MCD is 
not dramatically reduced when the field is decreased : 
this is in good agreement if compared to the observation 
of Wolos et aliii, but our observations (although with a 
low signal-to-noise ratio) tend to indicate a decrease and 
a change in shape of the MCD at small field. Note that 
the MCD is small indeed, and as such, may be sensitive 
to fluctuations of the trigonal crystal field, which can be 
different in a bulk crystal and an epitaxial layer. 



Cubic cristal field 
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1407 meV 


Trigonal cristal field 


SliVtri 


6.4 meV 


Spin-orbit 

coupling 

of 5 T 2 state 


Xeff — K\ 

P 

\d\ 


1.1 meV 
0.6 meV 
0.02 meV 


of 5 E state 


\d'\ 


0.015 meV 


Ham reduction factor 


K 


0.1 



TABLE I: Parameters of the effective Hamiltonian. 
In the perturbation treatment of the dynamic Jahn- 



Teller effect, the ratio A e ///A is equal to the Ham reduc- 
tion factor: hence k rs O.f . As K=e 3S / 2 , the Huang-Rhys 
factor is 5 =1.5, in reasonable agreement with the val- 
ues describing the intensity of the one-phonon replica in 
absorption and cathodoluminescence (0.7 and 1.9 respec- 
tively). 

The other parameters, p and d, govern two compo- 
nents of the effective Hamiltonian which are built through 
symmetry arguments. They can be estimated as second- 
order terms in the perturbation scheme. For instance, 
we expect a second-order contribution of the spin-orbit 
coupling through the excited vibrational states of the 
ground triplet, 29 , another contribution through the E 
doublet (the Xte term of Ref.l30l). which is opposite in 
sign, and a third contribution from spin-orbit coupling 
through excited triplet states and spin-spin coupling. 32 
Note however that this approach (perturbation treatment 
with coupling to a single E mode) is oversimplified and a 
full calculation, possibly involving also a coupling to J 1 2 
modest, would be necessary. This calculation, performed 
for (Zn,Cr)S and (Zn,Cr)Se, has even shown a stronger 
coupling to T2 mode than E mode. 

Nevertheless, the description of the magneto-optical 
data has been improved by introducing a simple descrip- 
tion of the dynamic Jahn- Teller effect. 

It may be also useful to define a spin Hamiltonian 
acting within the ground orbital singlet, which is sepa- 
rated by 6 meV from the doublet by the trigonal field. 
Its form for a spin 2 can be found in Ref. I'2("l3'jl3 ll In a 
wurtzite semiconductor, it remains : 



H = (9\\ 
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a(S 4 x 



102 N 
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As above, c labels the [000.1] axis of the wurtzite struc- 
tures, and (x, y, z) the cubic axes associated with the 
tetrahedron of nearest neighbors. 

The energy levels of the ground state calculated with 
the effective Hamiltonian (1) can be well reproduced us- 
ing the spin Hamiltonian (4), as shown by the red dot- 
ted lines in Fig. calculated for the parameters given 
in table II. We obtain g-factors slightly smaller than 2 
for both crystallographic directions (the orbital contri- 
bution is reduced by the Jahn- Teller effect), and a large 
axial anisotropy parameterized by D, which dominates 
the fine structure. As for Mn 24 "^ the cubic anisotropy 
terms a and F are much smaller than the axial term D. 
From the energy levels calculated in zero field with (1), a 
rough estimation of a — F w 0.03 meV can be obtained. 
The fine structure parameters have the same order of 
magnitude than the ones measured for Cr 24 ~(3d 4 ) ions in 
II-VFsi^ It should be noted that these parameters are 
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Mn 3+ (3d 4 ) S=2 



1.91 

1.98 

0.27 meV 

« 0.03 meV 



Mn 2+ (3d 5 ) S=5/2 



1.9994 

1.9994 
-2.5 fieV 
■ 0.07 neV 



TABLE II: Spin Hamiltonian parameters describing the 
ground spin multiplet of a Mn 3+ (3d 4 ) ion in GaN. The cor- 
responding parameters for a Mn 2+ (3d 5 ) in GaN— are given 
for comparison. 




1.4 1.5 1.6 1.7 
Energy [eV] 



12 3 4 5 6 
Mn composition [%] 



FIG. 7: (a) Absorption spectra for 0.06, 0.8, 1.7 and 5.7% Mn 
in GaN layers (b) area of the absorption band integrated over 
the energy range shown in (a), as a function of the Mn content. 
The samples with 0.03, 0.06 and 1.7 % Mn which have been 
studied in more details are shown by closed symbols. 



much larger than the corresponding ones determined by 
EPR for Mn 2 + ions in GaN (3d 5 configuration) 3 

We can also calculate the magnetization of isolated Mn 
ions. The result is not much different from the aver- 
age magnetization calculated under the assumption of 
a static Jahn- Teller effecti^ 5 - As expected, it is strongly 
anisotropic, the c-axis being a hard axis. 

The very sharp spectra described above have been ob- 
tained on samples with a very low Mn content, typically 
less than 0.1%. Very similar spectra, although broad- 
ened, have been observed up to 1% Mn (see an example in 
Fig- EH)- The intensity of the whole structure is propor- 
tional to the Mn content determined by SIMS (straight 
line in Fig. 0d). We could follow the absorption band 
in the near infrared up to 6% Mn, using FTIR in zero 
field. The strong broadening considerably smooths the 
structure, and finally it becomes difficult to eliminate the 
effect of the interferences and to determine the lineshape 
with a reasonable accuracy (top of Fig. EJu). Neverthe- 
less, the integrated area of the absorption band keeps 
increasing with the Mn content (Fig. [7Jd), showing that 
in most samples the majority of the Mn ions is in the d 4 
configuration. 



IV. MAGNETO-OPTICAL SPECTROSCOPY AT 
THE BANDGAP 

Magnetic circular dichroism at the band gap edge has 
been studied with the magnetic field applied parallel to 
the c-axis. For these studies, 0.5 /an thick (Ga,Mn)N 
epilayers have been grown on A1N MOCVD buffer lay- 
ers deposited on sapphire. A 100 nm thick GaN buffer 
layer had to be grown first in order to adjust the growth 
conditions. Here we describe the results obtained on a 
sample with 1.7% Mn. 

Fig.[H]shows the transmission for er + and o~ polarized 
light, measured at 1.7 K and 11 T using a linear polarizer 
and a Babinet compensator set as a quarter retarding 
plate. 

A broad absorption which rises with the Mn content is 
observed above 2.1 eVA It was attributed to a band- 
level transitionjSiiS At higher energy, the band edge rises 
in two steps: a weakly polarized one at low energy, which 
we shall ascribe to the (albeit thin) GaN buffer layer; a 
second one, with a clear circular polarization, 36 which 
we shall show to be due to the "giant Zeeman effect" 
in (Ga,Mn)N. Note the unusual polarization which is ob- 
served, with the a~ edge shifted to lower energy, in agree- 
ment with Rcf. HI 
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FIG. 8: (color online) Transmission of right- and left-polarized 
light (red dash line and blue solid line respectively) of the 
(Ga,Mn)N layer with 1.7% Mn, and magnetic circular dichro- 
ism spectrum at 11 T in the Faraday configuration (green 
solid line). The inset shows the MCD at 5 T with a broader 
energy scan. 



Fig.Elshows also the MCD, [I{a+) - I(cj-)]/[I(a + ) + 
I(a~)). The broad absorption band observed^ 6 , between 
2.1 and 3.5 eV shows no measurable polarization (inset 
of Fig.lHJ: MCD is only observed in resonance with the 
bandgap edge. 

If we assume that the absorption edge rigidly shifts in 
opposite directions in <r + and <j~ polarizations, the MCD 
can be expressed ast^i 




MCD 



A\n{T(E)) AE 
&E 2~' 



(5) 



3.40 3.45 3.50 3.55 3.60 
Energy [eV] 

FIG. 9: Magnetic circular dichroism divided by the effec- 
tive Zeeman splitting AE e ff, for the (Ga,Mn)N sample with 
1.7% Mn at 1.7 K. 



with T(E) the zero-field transmission measured with 
non-polarized light at energy E, and AE = E(a + ) — 
E(a~) the Zeeman splitting. In the absence of any iden- 
tified excitonic structure, we shall use this expression to 
extract an effective Zeeman splitting AE e ff. A minimal 
requirement is that the ratio MCD/AE e ff should not 
depend on the applied field. This is checked in Fig. [5] 
for the high-energy part of the MCD (3.50 to 3.55 eV) 
observed at different values of the applied field, which 
indeed exhibits a constant shape, so that we can extract 
a value of AE e ff(H). This value is plotted in Fig. ^3 
together with the value measured on a similar layer of 
pure GaN. In GaN, AE e ff(H) is linear in field and does 
not depend on the temperature. In (Ga,Mn)N, it in- 
creases non-linearly with the applied field, and decreases 
if the temperature increases. Above 100 K, it assumes 
the same value as in pure GaN. Note also in Fig. [5] that 
the low energy part of the MCD spectrum, at 3.49 eV, 
keeps increasing at high field: this part is actually char- 
acterized by a different effective Zeeman splitting (not 
shown), equal to that of GaN. Hence we attribute this 
signal to the thin GaN buffer layer. 

If the effective Zeeman splitting we observe contains a 
contribution due to the giant Zeeman effect, i.e., to the 
interaction between the bands of the semiconductor and 
the localized spins of Mn, it is expected to be related to 
the magnetization of the Mn system, and to exhibit a 
specific dependence on the Mn content and on the field 
and temperature. This we check now. 

Fig. Oik shows MCD spectra at 1.7 K and 11 T for 
(Ga,Mn)N layers with different values of the Mn content, 
from to 5%. The MCD peak shifts to the blue upon 
increasing the Mn content. It shows a strong resonance 
which markedly increases up to 2%, then one observes 
a broadening of the band edge and of the MCD signal, 
and a decrease of the MCD peak intensity. However, 
the effective Zeeman splitting AE e ff keeps increasing, as 
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FIG. 10: Magnetic field dependence at 1.7 K of the effec- 
tive Zeeman splitting AE e ff observed at the bandgap of pure 
GaN (open symbols) and of (Ga,Mn)N with 1.7% Mn (closed 
symbol) . 
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FIG. 11: (color online) (a) MCD spectra at 1.7 K and 11 T 
for (Ga,Mn) layers with different values of the Mn content 
from to 5%, the arrows indicate the photon energies used 
in Fig. EI]; (b, symbols, left scale) AE e ff(x) — AE e ff(0) as a 
function of the Mn content x ; (right scale) density of Mn ions 
(x, blue solid line) and of isolated Mn ions (x e ff, green line). 



shown in Fig. Illb . which displays the change of AE e ff 
due to the incorporation of Mn, measured at 1.7 K and 
11 T. The dependence on x is sublinear, and actually 
quite close to the variation of the density of free spins, 
x e ff, which would be expected if nearest-neighbor Mn- 
Mn pairs are blocked by a strong antiferromagnetic in- 
teraction. We shall come back to this point later. 

Fig. ^3 suggests a paramagnetic-like behavior of the 
Mn-induced change of AE e ff for the 1.7% Mn layer. For 
this sample, the giant Zeeman splitting can be deter- 
mined by subtracting the Zeeman splitting of the refer- 
ence GaN layer. Fig. ^| shows the giant Zeeman split- 
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FIG. 12: Giant Zeeman splitting AE eff (x, H, T)-AE eff (x = 
0,H), plotted versus the magnetic moment measured by 
SQUID between and 5T with the field applied along the 
c-axis (the line is a guide for eyes). The giant Zeeman split- 
ting is found to be proportional to the magnetization in the 
whole temperature and magnetic field range. 



ting, plotted versus the magnetic moment, measured by 
SQUID with the field applied along the c-axis. In the 
whole temperature (2K-20K) and field range (0-5T), the 
giant Zeeman splitting is found to be proportional to 
the magnetization. This, and the dependence on the 
Mn content, clearly evidence a coupling between the s,p- 
electrons of the host semiconductor and the (i-electrons 
localized on the Mn 3+ ions. 

Of course, a determination of the strength of the spin- 
carrier coupling would be highly valuable. This is a 
rather simple question in the case of Mn-based II- VI 
DMS's, involving spins of high symmetry (Mn 2+ ) in a 
zinc-blende semiconductor with moderate excitonic ef- 
fects. In the case of a d A magnetic impurity, the p — d 
coupling involves not only the usual spin-spin term, but 
also an orbital coupling which has been particularly stud- 
ied for Cr 2+ in ZnSei^&iS We may expect this term to 
be quenched by the Jahn- Teller factor. The influence of 
the excitonic structure however cannot be neglected, as 
shown by the analysis of the giant Zeeman effect mea- 
sured in (Zn,Co)Or& another wide-gap DMS with the 
wurtzite structure and a small spin-orbit coupling, or by 
the analysis of stress effects in GaN epilayers4^ In both 
cases, the electron-hole exchange interaction leads to an- 
ticrossings of the various excitons (A, B and C) of the 
wurtzite structure, so that the shift of the exciton can be 
quite different from the shift of the involved carriers. 

Because of the broadening of the excitons, and of the 
overlap with the contributions from the GaN buffer layer, 
excitonic transitions are not observed in the present lay- 
ers, neither in reflectivity nor in transmission. Hence a 
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quantitative analysis is out of reach. Nevertheless, the 
MCD observed in transmission can be tentatively at- 
tributed to the exciton of lowest energy, which in GaN is 
the so-called A exciton of the wurtzite structured The 
giant Zeeman splitting of the corresponding electron-hole 
pair, in the Faraday configuration with the field along the 
c-axis, is^fi 

AE = E((T+) - E(<j-) =N (a-f3)x <S Z >, (6) 

where a and [3 are the s- and p-d coupling constants, 
respectively, N x is the spin density, and < S c > is the 
spin component of Mn along the c-axis (=-2 at satura- 
tion). From the spectra shown in Fig. [S] a negative 
sign of a — f3 is inferred. This is opposite to what was 
found in most DMS's studied so far, but in Cr-doped 
II- VI semiconductors^^ another realization of the d 4 
configuration. As a is usually positive, and small with 
respect to /?, that means a ferromagnetic interaction be- 
tween the holes in the valence band and the localized 
spins, as it has been predicted for DMS's with a less 
than half filled d-shell>22 

Finally, it should be stressed that, in the absence of a 
complete identification of the various excitons, a quan- 
titative analysis of the giant Zeeman splitting would be 
adventurous. In particular, we expect the excitonic shifts 
to show a complex dependence on the magnetization, 
with anticrossings as soon as the exchange splitting of 
the bands becomes of the order of the zero-field exci- 
tonic splittingjiMi That will cause the MCD to satu- 
rate, which could contribute to the saturation observed in 



Figlllb as a function of the Mn content, or in Fig^|as a 
function of the applied field (while the magnetization cal- 
culated for a single Mn ion is still far from saturation); 3 - 



V. CONCLUSION 

A detailed study, by magneto-optical spectroscopy in 
the near infrared, of wurtzite (Ga,Mn)N layers grown 
by plasma-assisted molecular beam epitaxy, confirms the 
substitution of Ga by Mn with the d A configuration 
(Mn 3+ ). Very sharp spectra obtained on samples with 
a low Mn content, and cathodoluminescence spectra, are 
reasonably well understood using a crystal field model 
with a dynamic (not static) Jahn- Teller effect. 

Samples with a slightly larger Mn content show a MCD 
signal in transmission, resonant with the bandgap edge. 
The effect increases with the Mn content and behaves as 
the measured magnetization. This evidences a coupling 
between the Mn 3+ rf-electrons and the carriers of the host 
semiconductor (" giant Zeeman effect" ) ; the sign of the 
MCD agrees with a ferromagnetic p-d interaction. 
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